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The Stirling numbers of the first kind can be represented in terms of a new class of 
polynomials that are closely related to the Bernoulli polynomials. Recursion relations 
for these polynomials are given. 

The Stirling numbers have wide application in science and engineering. For example, the 
generating function of the Stirling numbers of the first kind appears in the inversion of a 
Vandermonde matrix Q , which is used in curve fitting, coding theory, and signal processing 



In this paper we show that there is a connection between the Stirling numbers of the 
first kind and the Bernoulli polynomials. The unsigned Stirling numbers of the first kind 
(— l) n ~ m Sn are defined as the number of permutations of n symbols which have exactly m 
permutation cycles 0,0]. Let us construct the positive integers T n fc by 
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Here is a table of the numbers T n 
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Note that the horizontal sums of the numbers in this table are factorials. That is, 1 + 1 = 2!, 
1 + 3 + 2 = 3!, 1 + 6 + 11 + 6 = 4!, and so on. 

Let us construct a new set of polynomials Pk{n) that reproduces these Stirling numbers. 
For example, the numbers (0, 1, 3, 6, 10, . . .) in the first column in this array as a function 
of the row label n are given by the polynomial 



where 
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Pn(n) = 1. 
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The numbers (0, 0, 2, 11, 35, . . .) in the next column in this array as a function of n are 
given by 

_L n(n _l) (n _ 2)P 1 (n), 
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where 



Pi(n) 
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The numbers in the next two columns in this array are given by 
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In general, the numbers along the fcth column are expressible in the form 
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The next seven polynomials Pk{n) are 
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The structure of these monic polynomials strongly resembles that of the Bernoulli poly- 
nomials B m (x) First, like the Bernoulli polynomials, every other polynomial in (j2J) - © 
has a constant term. Second, recall that the constant terms in the Bernoulli polynomials are 
Bernoulli numbers, i?2m(0) = -82m, and observe that the constant terms in the polynomials 
-p2m-i(^) are expressible in terms of Bernoulli numbers: 
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Furthermore, both the polynomials P m (n) and the Bernoulli polynomials B m (x) satisfy 
very similar recursion relations. The polynomials P m (n) satisfy 



P 2m (n) = (2m + l)nPam-i(ra) + X^ -2 ) 
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while the Bernoulli polynomials obey the same recursion relation with an inhomogeneous 
term: 

B 2m+l {x) = (2m + l)xB 2m {x) + V(-2) fc - - ",\ ) , x fc+1 B 2m _ k (x) 
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Finally, we note that P m (n) satisfies a recursion relation in terms of the coefficients of 
the Bernoulli polynomials: 
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P m+ i{n) = nP m (n) ^— V 4 fe P m +i_ 2 fc(^) [coefficient of x m+2 2fc in S m+2 (x)] 
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